Published online 19 April 2011 



Nucleic Acids Research, 2011, Vol. 39, No. 14 5991-6001 

doi:10.1093/nar/gkr219 



Evidence for an evolutionary antagonism between 
IVIrr and Type III modification systems 

Mehari Tesfazgi Mebrhatu\ Ewa Wywial^, Anirban Ghosh\ Chris W. Michiels\ 

Ariel B. Lindner^, Fran9ois Taddei^, Janusz M. Buinicki^'"*, Laurence Van Melderen^ and 

Abram Aertsen^'^'^'* 

^Laboratory of Food Microbiology and Leuven Food Science and Nutrition Researcin Centre (LFoRCe), 
Centre for Food and IVIicrobial Teclnnology, Department of IVIicrobial and IVIolecular Systems (M^S), Faculty of 
Bioscience Engineering, Katholieke Universiteit Leuven, Kasteelpark Arenberg 22, B-3001 Leuven, Belgium, 
^Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in 
Warsaw, ul. Ks Trojdena 4, 02-109 Warsaw, Poland, ^Laboratoire Genetique Moleculaire Evolutive et Medicale, 
INSERM U571, Universite Rene Descartes-Paris V, rue du Faubourg St. Jacques 24, 75014 Paris, France, 
"^Bioinformatics Laboratory, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam 
Mickiewicz University, ul. Umultowska 89, 61-614 Poznan, Poland and ^Laboratoire de Genetique et de 
Bacterienne, Institut de Biologie et de Medecine Moleculaires, Universite Libre de Bruxelles, rue des 
Professeurs Jeener et Brachet 12, 6041 Gosselies, Belgium 

Received March 9, 2011; Revised IVIarcPi 23, 2011; Accepted IVlarcli 24, 2011 



ABSTRACT 

The Mrr protein of Escherichia coli is a laterally 
acquired Type IV restriction endonuclease with spe- 
cificity for methylated DNA. While Mrr nuclease 
activity can be elicited by high-pressure stress in 
E. coli MG1655, its (over)expression per se does 
not confer any obvious toxicity. In this study, 
however, we discovered that Mrr of E. coli MG1655 
causes distinct genotoxicity when expressed in 
Salmonella typhimurium LT2. Genetic screening 
enabled us to contribute this toxicity entirely to the 
presence of the endogenous Type III restriction 
modification system (StyLTI) of S. typhimurium 
LT2. The StyLTI system consists of the Mod DNA 
methyltransferase and the Res restriction endo- 
nuclease, and we revealed that expression of the 
LT2 mod gene was sufficient to trigger Mrr activity 
in E. coli MG1655. Moreover, we could demonstrate 
that horizontal acquisition of the MG1655 mrr locus 
can drive the loss of endogenous Mod functionality 
present in S. typhimurium LT2 and E. coli ED1a, and 
observed a strong anti-correlation between close 
homologues of MG1655 mrr and LT2 mod in the 
genome database. This apparent evolutionary an- 
tagonism is further discussed in the light of a 
possible role for Mrr as defense mechanism against 
the establishment of epigenetic regulation by foreign 
DNA methyltransferases. 



INTRODUCTION 

In general, DNA restriction endonucleases (REases) team 
up with DNA methyltransferases (MTases) to form dedi- 
cated restriction modification (RM) systems, in which the 
MTase modifies DNA at specific sequences in order to 
protect it from cleavage by the cognate REase (1,2). 
Because of this interplay, bacterial RM systems behave 
as a primitive immune system in which the REase restricts 
the cellular entry of DNA lacking the proper methylation 
signature. As such, it can protect the host against deleteri- 
ous phage infection or lateral acquisition of foreign DNA 
(3). However, as the genes encoding such RM systems are 
typically closely linked and often organized into an 
operon, they can also behave as selfish addiction modules 
that counteract their loss from the cell by post- 
segregational kilhng (4,5). In this phenomenon, sibhngs 
that lost the RM module will experience a decline in cor- 
responding MTase activity, and eventually fail to protect 
all genomic target sites from cleavage by the residual 
cognate REase. 

Depending on their operational characteristics, RM 
systems can be categorized into three main types (Type 
I-IIl) (6). As such. Type I systems are unique in requiring 
a dedicated specificity protein that directs activity of the 
corresponding MTase and REase to the target sequence, 
with cleavage itself occurring at variable distances from 
this site. In contrast, Type II MTase and REase pairs in- 
dependently recognize the same target sequence and most 
often function as separate entities, with cleavage occurring 
at a defined position within or close to the recognition site. 
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Finally, Type III systems display specificity for a short 
asymmetric recognition sequence, and cleave at a charac- 
teristic distance away from it. Importantly, binding to the 
target site is strictly mediated by the MTase, so that 
restriction is performed by a MTase-REase complex 
rather than by the REase alone. While all RM systems 
require AdoMet for donation of the methyl group. Type 
I and III systems additionally depend on ATP for their 
activity (6). 

Interestingly, however, a small number of REases 
operate in the absence of a cognate MTase and have 
been grouped into a fourth category (Type IV). These 
REases invariantly bear a distinct specificity for modified 
(e.g. methylated, hydroxymethylated or glucosyl- 
hydroxymethylated) DNA, although their recognition 
site is often poorly defined (6). Well-known examples 
of such sohtary and modification-dependent REases in 
Escherichia coli are McrA, McrBC and Mrr, of which 
the corresponding genes are typically located on laterally 
acquired genetic elements such as el4 (mcrA) and the 
'Immigration Control Region' (mcrBC and mrr) (7-9). 

While McrA and McrBC are relatively well- 
characterized, and shown to be involved in restricting 
phage infection (10-14), the physiological role of Mrr in 
the cell has been less obvious so far. Recently, however, 
our interest in Mrr was raised as it appears to be a driving 
part of the peculiar piezophysiology of E. coli. More 
specifically, it was shown that cellular filamentation and 
prophage induction caused by sub-lethal high-pressure 
shock in E. coli K12 were the result of a RecBCD- 
dependent SOS response that was instigated by Mrr endo- 
nuclease activity (15-18). Unfortunately, it so far remains 
unclear how exactly Mrr becomes activated by a physical 
stress such as high pressure. Further study, however, 
revealed that not all Mrr homologues behave the same 
way, and the Mrr protein encoded by Salmonella 
typhimurium LT2, in fact, could not be activated by 
high-pressure shock (19). 

The current study starts with the serendipitous observa- 
tion that the Mrr protein of E. coli MG1655 causes 
distinct genotoxicity when expressed in 5. typhimurium 
LT2, and eventually results in postulating an evolutionary 
antagonism between Mrr and Type III MTases. 

MATERIALS AND METHODS 

Strains and construction of mutants 

Salmonella typhimurium LT2 (20, laboratory collection), 
Escherichia coli KI2 MG1655 (21, laboratory collection) 
and E. coli EDI A (22, kindly provided by Eric Denamur, 
Institut Pasteur, Paris, France) were used as parental 
strains in this study. 

An LT2 derivative cured of its four active prophages 
(i.e. strain MA8508 constructed by Lemire and Bossi, 
manuscript submitted for publication) was kindly provided 
by Lionello Bossi (Centre National de la Recherche 
Scientifique, Gif-sur-Yvette, France), and is referred to 
as LT2Acp in this study. Strain LT2K2 is a derivative 
of LT2 and harbors a MudK (23) based translational 
lacZ fusion to the SOS-responsive STM0912 gene. 



and its use as an SOS indicator strain has been validated 
previously (19). LT2 recA was constructed by cotransduc- 
ing the recAl and «-/-202::7>j10 alleles from TT521 (kindly 
provided by John Roth, University of CaHfornia at 
Davis, USA) with P22 HT105/1 int-201 (24). In strain 
LT2 AStyLTI, the StyLTI mod res operon was deleted 
based on the protocol of Datsenko and Wanner (25), 
and using an amphcon prepared on pKD4 with the fol- 
lowing primers: 5'-ATAAGAAATTGATTGCTGCACC 
TTAAAATTATCAACGGAGCGAAGATATCGTGTA 
GGCTGGAGCTGCTTC-3' and 5'-ATTTTATATGGTC 
TATCAAAAAGATAACCGTTAAGGTTCATTCAAA 
CCCGCATATGAATATCCTCCTTAG-3'. This proced- 
ure replaced the mod res operon with an /;-? flanked kana- 
mycin resistance cassette, which could subsequently be 
excised by transiently equipping this strain with plasmid 
pCP20 expressing the Flp site-specific recombinase (26). 
Strain LT2K2 Amod res was subsequently constructed 
by transducing STM0912::MudK to LT2 Anwd res 
using P22 HT105/1 int-201. 

Strain MG1655 mrr harbors an mrr::Tn5 allele, which 
was obtained through an in vitro transposomics procedure 
using the EZ-Tn5 transposome kit (Epicentre, Landgraaf, 
The Netherlands) and an amplicon of the MG1655 mrr 
gene (mrr^°^^^^). The latter was obtained by PCR using 
the foUowing primers: 5'-CGATTCTAGACTGCAAGG 
ATGTACTATGACG-3' and 5'-CGATAAGCTTGCGT 
TTGCGGGGTTGAGG-3'. The resulting mrrv.TnS allele 
confers kanamycin resistance and was exchanged with the 
wild-type mrr allele according to the protocol of Datsenko 
and Wanner (25). Strain MG1655 XsfiAv.lacZ was con- 
structed by lysogenizing strain DJ624 (MG1655 
AlacX14 mat.-.lacf-) (27) with XsfiAv.lacZ, originating 
from CSH50 ksfiAwlacZ (28). Subsequently, Strain 
MG1655 ksfiAwlacZ mrr was constructed by transducing 
the mrr.:Tn5 allele from MG1655 mrr into MG1655 
XsfiA--lacZ with Plv/r (29). 

Strains were transformed with the necessary plasmids 
by electroporation, while curing of temperature-sensitive 
plasmids, such as pKD46 and pCP20, was performed by 
growing the corresponding strain at the non-permissive 
temperature in the absence of plasmid selection and sub- 
sequently isolating a clone that had lost the plasmid. 

Growth conditions 

Stationary phase cultures were obtained by growth in 
lysogeny broth (LB) (30) for 21 h at 37°C under 
wefl-aerated conditions. Late exponential phase cultures 
in turn were prepared by diluting stationary phase 
cultures 1/100 or 1/1000 in fresh pre-warmed LB, and 
allowing further incubation at 37°C until OD^oo = 0.6 as 
described earlier (31). When appropriate, the following 
chemicals (Applichem, Darmstadt, Germany) were 
added to the growth medium at the indicated final con- 
centrations: ampicilhn (100|ig/ml), chloramphenicol 
(30 ng/ml), kanamycin (50 |ig/ml), glucose (0.2 or 0.02 
%) and arabinose (0.2 or 0.02%). 

Where indicated, quantitative measurement of growth 
was performed with either a Multiskan Ascent reader 
(Thermo, Aalst, Belgium) or a Bioscreen C Microbiology 
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Reader (Labsystems Oy, Helsinki, Finland), which automat- 
ically measure OD^oo of cultures growing in thermostatted 
microplate wells at regular time intervals. 

Construction of plasmids 

Plasmid pFPVlS-P^^^j-morf'^^^ is derived from the 
pFPV25 backbone (32), confers ampicilhn resistance and 
conditionally expresses Mod^^^ in the presence of arabin- 
ose. For its construction, an amplicon of the mocf'^^ allele 
was first obtained by PCR on genomic DNA of LT2, 
using the following primers: 5'- AAGTTCTAGAGGAG 
CGAAGATATCATGTTG -3' and 5'- CGTAAAGCTTA 
TCATCCCCTCTTAATAACG -3'. Subsequently, this 
amplicon was digested with Xbal and Hindlll, and 
cloned in ^¥¥V 25-? sAD-gfp (i-e. pAAlOO; 19) digested 
with the same enzymes to replace the gf'p gene and yield 
pFYVlS-F BAD-mocF^"^ . Similarly, plasmid pFPV25-Pb^£,- 
mocF^^^ was constructed by preparing a PCR ampHcon 
of the mocF^^^ allele from genomic DNA of EDI A, using 
the following primers: 5'- ACGC TCTAGA GATACAGA 
TACGTGGAA-3' and 5'- CGTA CTGCAG CATCTCAT 
TCCCTTACAG-3'. Subsequently, this amplicon was 
digested with Xbal and PstI, and cloned in pFPV25- 
PBAD-gfp digested with the same enzymes to replace the 
gfp gene and yield nFPV25-P BAD-mocf-°^^. Plasmid 
pFPW25-P BAD-mrr'^° was constructed similarly, but 
using the OTn-MG1655 amplicon obtained as described 
above. Plasmid pFPV25-P5^£)-7iirr^^^ was constructed 
earlier in a similar way (19). Where necessary, plasmid 
pFPV25-P£^£, served as a control, and was obtained by 
removing the gfp gene from pAAlOO with Xbal and 
Hindlll, and reclosing the resulting fragment by 
self-ligation after blunting. 

Plasmid pACYCm-P BAD-mrr^°^^^^ is derived from 
the pACYC184 backbone (33), confers chloramphenicol 
resistance and conditionally expresses Mrr^°"' ^ in the 
presence of arabinose. For its construction, a fragment 
encoding P^^/j-wirr'^^'*^^ was cut out from the 
pFPV25-P BAD-mrr'^°^^^^ vector with EcoRI and 
Hindlll, and subsequently blunted. This fragment was 
then ligated into the pACYC184 vector, which was first 
opened with Hindlll and BamHI and blunted, to yield 
pACYC184-P5^B-'«'-'-'^°'"^ Plasmid pAA810, which is 
also pACYC184 based, was constructed earlier (18) and 
harbors the «in-'^°"'^^ locus under the control of its own 
promoter. 

Screening for suppression of Mrr'^^'*^^ lethality in 
LT2 recA 

After we established that pACYClM-P BAD-mrr^°^^^^ 
conferred severe lethality in LT2 recA when plated on 
M9 minimal medium (30) containing 0.2 % arabinose, 
two different screens were initiated to find suppressors 
of this lethahty. In a first screen, LT2 recA pACYC184- 
PsAD-'Tirr^'^^^ ^ was transformed with a hbrary contain- 
ing random fragments of MG1655 genomic DNA cloned 
in pBR322 (kindly obtained from Susan Gottesman, 
National Institutes of Health, Bethesda, USA). In a 
second screen, a random Mud} knock-out hbrary was 
constructed in LT2 recA pACYCm-P BAD-mrr'^'^^^^^ 



using the MudJ transposon delivery system described by 
Hughes and Roth (23). Next, the corresponding plasmid 
and transposon libraries were plated on M9 medium com- 
plemented with 0.2% arabinose and appropriate antibiot- 
ics, and clones able to grow on such plates were isolated 
and purified. Subsequently, the corresponding pBR322- 
plasmids were transformed to a fresh LT2 recA 
pACYC184-P5^£,-nwT^^"'^^ strain, while the correspond- 
ing transposon insertions were transduced to a fresh LT2 
pACYC184-Ps^£,-ww7-^^^'''-'^' strain after which the recAl 
allele was co-transduced into the resulting strain as 
described above. These resulting strains were examined 
to check whether the Mrr'^*^'*^'^ resistant phenotype 
could be confirmed. 

One valid Mud] clone arose from this selection proced- 
ure, and its insertion site was determined by subcloning 
the chromosome of the corresponding mutant in pUC18 

(34) and selecting for the kanamycin resistance marker 
encoded on the Mud] transposon. 

Measuring of SOS induction and DNA damage 

SOS induction was measured via standard Miller assays 

(35) , using LT2K2 or MG1655 XsftA::lacZ-hased indicator 
strains. For the biochemical detection of DNA breaks and 
their distribution throughout a population, terminal 
deoxyribonucleotide transferase dUTP nick end labehng 
(i.e. TUNEL assay; in situ CeU Death Detection Kit, 
Fluorescein, Roche Mannheim, Germany) was used, 
based on the manufacturer's instructions and the amend- 
ments published earlier (36,37). Briefly, 1 ml of a culture 
was washed two times with phosphate-buffered saline 
(PBS; 0.14M NaCl, 2.7 mM KCl, 1.5 mM KH2PO4, 
9.6 mM Na2HP04, pH 7.3), and fixed with ice-cold 4% 
paraformaldehyde for 30min at room temperature. Cells 
were then pelleted and resuspended in ice-cold permeabil- 
ization solution (0.1% Triton X-100, 0.1% sodium citrate) 
and incubated for 2min at 4°C. Again ceUs were washed 
one time with PBS and labeled using TUNEL reagents 
for 1 h in the dark at 37°C. Finally, ceUs were washed 
two times with PBS and analyzed using flow cytometry 
(Guava Easy Cyte Plus, Guava Tehnologies, Inc. 
Hayward, CA, USA). To generate a positive control, an 
aliquot of fixed and permeabilized control cells was first 
exposed to the blunt 4-cutter BsuRl for 1 h at 37°C, before 
being washed and labeled with TUNEL reagents. 

Constructing the phylogenetic tree of E. coli and 
S. enterica strains 

To construct the phylogenetic tree of E. coli and 
S. enterica strains, all full-length Escherichia and 
Salmonella 16 S rRNA gene sequences were downloaded 
from the Greengenes database core set (38) in the FASTA 
aligned format. Next, E. coli EDI A 16S rRNA gene 
sequence was identified using NCBl nblast (with E. coli 
CFT073 16S rRNA as a seed) and added manually to the 
dataset. Gene sequences that could not be successfully 
aligned originated from partially sequenced genomes or 
were duplicative, and were removed from the dataset. 
Next, the dataset was re-aligned using Greengenes (38), 
which ahgns 16S rRNA gene sequences to 7682 characters 
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full-length gene templates. All thresholds were kept at 
default values. Finally, based on the ahgnment, the phylo- 
genetic tree was calculated using MEGA 4.0 (39) employ- 
ing the Minimum Evolution method and assuming a 
Jukes-Cantor model of nucleotide substitution. 
Bootstrap values based on 1000 replications are hsted as 
percentages at the branching points. 

Calculating the anti-correlation of Mrr"^^'^^^ and 
Mod^^^ homologues 

Sequences of E. coli Mrr (gi: 127320) and 5. typhimurium 
Mod (gi: 300193) were used as queries for PSI-BLAST 
(40) searches (£- value < le~^) against the nr database 
using default parameters and run until convergence. We 
retrieved 2001 Mrr and 6878 Mod sequences. After 
removing sequences from partially sequenced genomes 
and redundant sequences, we performed sequence cluster- 
ing based on pair-wise BLAST similarity scores, using 
Cluster Analysis of Sequences (CLANS) (41). The cluster- 
ing was completed at /" = 0.012 for Mrr and P = 0.004 for 
Mod, respectively. The reported /"-values give a 
weU-resolved separation of multiple distinct clusters. The 
/•-value was chosen empirically, given the /'-value plot for 
each data set, which shows a histogram of the number of 
sequences for each is-value below a given value. For 
example, a cut-off of 0.004 will exclude connections 
worse than 0.004. Clusters are thought of as robust if 
smaU changes in cut-off values do not result in major 
changes in their content. Next, we extracted members of 
clusters containing the Mrr query sequence (gi: 127320) 
and the Mod query sequence (gi: 300193) for correlation 
analysis. At the end, we obtained 75 Mrr'^°'*" and 211 
Mod^^"^ homologues derived from 272 fully sequenced 
genomes and eliminated duplicate species to obtain 45 
]^j.j.MGi655 J 5^ Mod^'^ homologues derived from 
192 fully sequenced genomes. 

The Pearson r correlation coefficient measures the 
degree to which values of two variables are linearly 
related to each other. It is defined as the covariance of 
two variables divided by the product of their standard 
deviations and was calculated for the Mrr'^°"'^^ and 
Mod^'^^ families using STATISTICA 8 (StatSoft, Inc., 
Tulsa, OK, USA). We also determined the probabihty 
that the observed correlation is real and not a chance 
occurrence. The obtained correlation of —0.52 for 45 
]y[j.j.MGi655 ^j^j Mod^^^ homologues, derived from 
192 fully sequenced genomes, is less than the critical 
value for df= 190, a = 0.05, two-tailed test with 
P< 0.0001. 

Determination of the horizontal transfer index for 
homologs of mrr'^^'^^^ and mod''^^ 

Nakamura et al. (42) have pubhshed an analysis of 116 
completely sequenced prokaryotic genomes, in which they 
calculated an index of potential distant horizontal transfer 
for all genes, by comparing the frequency of 'words' of 
pentanucleotide length within each gene with the average 
word frequency of the entire genome. This index can be 
used to detect only recent transfers of genes from species 
with distinct gene sequence compositions, because genes 



transferred between similar genomes usually display 
similar sequence compositions and genes transferred a 
long time ago tend to accommodate their sequence com- 
position to that of the new host. We examined the hori- 
zontal transfer index of genes encoding Mrr"^*^'^^^ and 
Mod^^^ homologues in the updated set of 165 genomes 
(Yoji Nakamura and Takashi Gojobori, personal 
communication) . 

RESULTS 

Expression of Mrr"^^'^^^ is harmless in E. coli MG1655 
but genotoxic in S. typhimurium LT2 

Previously, we have shown that, in contrast to MG1655, 
LT2 is unable to mount the SOS response after high- 
pressure shock (19). However, during our initial attempts 
to possibly restore this feature by complementing LT2 
with a plasmid encoding the under the 

control of its own promoter (i.e. pAASlO), we obtained 
colonies with a sick translucent appearance. Moreover, 
when transformants were grown in liquid broth, lysis of 
the culture (Figure lA) and concomitant accumulation of 
cellular debris could be clearly observed. 

As the genome of LT2 contains four functional pro- 
phages (20), we postulated that cell lysis was caused by 
prophage activation due to Mrr'^°'*'^^-mediated DNA 
damage. Upon further investigation, an LT2 mutant 
deprived of its active prophages (LT2A(p) indeed failed 
to show lysis (Figure IB). Moreover, we could clearly 
detect SOS induction (Figure 1 C) and even confirm phys- 
ical DNA damage (Figure ID) upon expression of 
]yjj.j.MGi655 LJ2, both rcvcaling the severe genotoxicity 
of this protein in LT2. 

Despite its constitutive toxicity in LT2, it is important 
to note that Mrr'^*^'*^^ was harmless upon similar ectopic 
expression in its natural MG1655 host (data not shown), 
unless such strains were exposed to high-pressure shock 
(18). We, therefore, assumed this surprising dissimilarity 
to stem from differences in the physiology of both strains, 
resulting either in dedicated mechanisms of Mrr'^'^'^^^-re- 
striction alleviation in MG1655 or dedicated mechanisms 
of Mrr"^*^'"-^ activation in LT2. 

The endogenous StyLTI Type III restriction modification 
system is responsible for Mrr'^^*^^^ toxicity in 
S. typhimurium LT2 

In order to find out why Mrr'^*^"'^^^ behaves so differently 
in closely related species, a selection scheme was devised 
in which we tried to suppress the lethality of Mrr'^°"'^^ in 
LT2. First of all, delivery of Mrr"^*^'*^^ was made depend- 
ent on the presence of arabinose (using pACYC184-P5^£)- 
^,.^MGi655-j^ in order to clearly define the onset 
of genotoxicity and avoid trivial suppressor mutations in 
the mrr^'^"'^ locus. As such, conditional expression of 
]yjj.j.MGi655 j.g(jy^,g(j viability of LT2 wild-type with 
almost 99.99% (i.e. ca. 4-log reduction), although 
this lethahty could be mitigated in the LT2A(p strain 
(Figure 2). Subsequently, in order to avoid picking up 
suppressor mutants that would trivially be affected in 
prophage induction, we used an LT2 recA mutant as the 
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Figure 1. Mrr"'^"'^^-mediated genotoxicity in S. typhimurium LT2. Growth of LT2 wild-type (A) or LT2A(p (B) equipped with either pACYC184 
(mrr—) or pAASlO (mrr+). (C) Activation of the SOS response in late exponential phase cultures of LT2K2 equipped with either pACYC184 (mrr—) 
or pAASlO (mrr+). P-Galactosidase activity in the LT2K2 background serves as a measure of SOS induction, and is expressed as Miller Units. (D) 
Flow cytoinetry of TUNEL-labeled late exponential phase cultures of LT2 pACYC184-PB.jo- 

,„,,MG1655 i„ 

LB in the presence of 0.02% glucose 

(mrr—) or 0.02% arabinose (mrr+). Fluorescence is a measure of DNA damage, and permeabiUzed cells of LT2 treated with the blunt four-cutter 
i?.vwRI to deliberately introduce DNA breaks in vitro were used as a positive control. 
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Figure 2. Mrr"'^"'"-mediated lethality in S. typinmurium LT2. Survival of indicated LT2 mutants, carrying pACYC184-P5^o-;;!)7-"'^"''^ and grown 
to stationary phase in LB, when plated on M9 minimal medium agar plates supplemented with 0.2% glucose (mrr—) or 0.2% arabinose (mrr+). Cell 
counts were expressed as Log N, in which N represents cfu/ml. 



latter does not support prophage activation but proved to 
be similarly hypersensitive to Mrr'^°'^^^ because of its 
inability to repair DNA damage (Figure 2). Finally, sup- 
pression of Mrr'^'^'*-^^^ toxicity in LT2 recA pACYC184- 
^BAD-f'^rr^'^^^^^ was attempted by introducing (i) random 
transposon knock-outs or (ii) a multi-copy plasmid hbrary 
of the MG1655 genome. As such, we hoped to identify cel- 
lular functions that either trigger or mitigate Mrr"^*^'^^^ 
activity in LT2, respectively. 

Of both strategies, only the knock-out approach 
resulted in a bona fide Mrr'^°"'^^ resistant mutant, and 
subcloning enabled us to map its transposon insertion site 
to the LT2 res gene (Figure 2). This gene is located imme- 
diately downstream of the mod gene, and the mod res 
operon encodes the Type III StyLTl restriction- 



modification system of LT2. In order to unequivocally 
corroborate the interference of the StyLTl system with 
]yjj.j.MGi655^ clcau Amod, Ares and Amod res mutants 
were constructed de novo in wild-type LT2 and 
examined. Confirmingly, Mrr'^°"'^^ lethality (Figure 3A) 
and genotoxicity (Figure 3B) proved to be completely 
abolished in LT2 Amod and Amod res mutants, and to 
be greatly attenuated in an LT2 A res mutant. These 
results clearly delineate the Type III StyLTl system as 
the sole feature responsible for constitutive Mrr'^*^'^^^ 
toxicity in LT2. 



Functional incompatibility between Mod and Mrr 



4LT2 



Based on the above results, it was postulated that Mod 
MTase activity would render the chromosome susceptible 
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Figure 3. The S. typhimitrium StyLTI RM system is responsible for Mrr'^'^'^''^-mediated toxicity. (A) Survival of indicated LT2 mutants, carrying 
pACYC184-PBy(D-mn''^'''^^' and grown to stationary phase in LB, when plated on M9 minimal medium agar plates supplemented with 0.2% glucose 
(mrr—) or 0.2% arabinose (mrr+). Cell counts were expressed as Log N, in which N represents cfu/ml. (B) Activation of the SOS response in late 
exponential phase cultures of indicated mutants of LT2K2 equipped with either pACYC184 (mrr—) or pAASlO (mrr+). P-Galactosidase activity in 
the LT2K2 background serves as a measure of SOS induction, and is expressed as Miller Units. 



to degradation by Mrr"^*^'^^^, and the corresponding genes 
were conditionally co-expressed in MG1655 mrr. While 
expression of mocf'^^ or m;v'^°'*'^^ alone was completely 
harmless in this background, their co-expression resulted 
in poor growth (data not shown) and massive SOS induc- 
tion (Figure 4). These data confirm the apparent incompati- 
bihty between both proteins, and identified Mod^^^ 
activity as instigator of Mrr"^*^'*^^ restriction. 

In order to examine whether this functional incompati- 
bility could also extend into mutual exclusiveness between 
mrr^'^"'^^ and mod^^^ loci, LT2 was equipped with 
pAASlO and 10 independent clones were followed for a 
number of generations. After only two passages, cultures 
grew without lysis (Figure 5A) or the accumulation of cel- 
lular debris, and gave rise to healthy colonies upon plating. 
Importantly, four out of these 10 independent suppressor 
mutants carried mutations in mod^^^ (Figure 5B), while 
others were affected in the mrr^'^^^^^ locus encoded on 
pAASlO. Interestingly, two of the compromised mod'^^~ 
alleles suffered frame-shift mutations within short homo- 
polymeric tracts present in the open reading frame 
(Figure 5B). These data indicate that Mrr'^°'^^^^ and 
Mod^^^ activities are mutually exclusive, and that acqui- 
sition of one could lead to functional inactivation of the 
other. 
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pACYC-Pg^-mrr"°i'S55 + + 

pFPV25-PBAo-'no^^™ + + 

Figure 4. Expression of the Mod^^^ MTase triggers Mrr"*"'"'^ activity 
in E. coli MG1655. Activation of the SOS response in late exponential 
phase cultures of MG1655 mrr XsfiAwlacZ equipped with the indicated 
combination of plasmids and grown in the presence of 0.02% arabin- 
ose. p-galactosidase activity in the MG1655 XsfiAwlacZ back- 
ground serves as a measure of SOS induction, and is expressed as 
Miller Units. 
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1121 CGAATAAAGAGGGGAAAAAGG 1141 
GGAATAAAGA GGGAAAAAGG 



1296 TGAAGAAGATGGGGGGCAGCG 1316 
TGAAGAAGAT GGGGGCAGCG 



870 AGGGAAATGGGGACCGGTGCC 890 
AGGGAA CGGTGCC 

******** 



898 GTGTCAAAAGGGCTGAGCGGC 918 
GTGTCAAAAGAGCTGAGCGGC 



Time (h) 

Figure 5. Spontaneous alleviation of Mrr'^'^"^^-mediated toxicity. (A) Typical restoration of normal growth of LT2 pAASlO during consecutive 
passages (1/100 dilution at indicated time points) in LB medium. Representative data are shown. (B) Mutations identified in the mod^^^ alleles of 
four independent LT2 pAASlO mutants in which Mrr"^*"'""^^ toxicity was suppressed. Numbers indicate the position of the base at the beginning and 
end of the displayed wild-type (i.e. upper) sequence with respect to the start of the open reading frame. Asterisks highlight the actual deletions or 
substitution in each of the four alleles (i.e. lower sequence). 
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Figure 6. Minimal evolution phylogenetic tree of 16S rRNA derived from fully sequenced Escherichia coli and Sahnonella enterica genomes. The 
number at each node corresponds to the percent bootstrap value (1.000 bootstraps). Additionally, NCBI Taxonomy ID as well as Mrr and Mod 
protein content is indicated in the lineage of each strain. The branches are colored blue for genomes with Mrr homologues, red with Mod homo- 
logues, green with both Mrr/Mod homologues, and yellow with no Mrr/Mod homologues. 



In the same context, it was noteworthy to find that close 
homologues of mrr^'^^^^^ and mocl^ ^ often seem to 
co-exist within S. enterica (Figure 6), which could impli- 
cate the functional loss of one of these alleles. Indeed, 
most S. enterica mrr alleles very closely resemble that of 
strain LT2, and we were unable to observe any Mod"^^^- 
dependent (or in fact any) genotoxicity upon deliberate 



over-expression of Mrr in LT2 (data not shown), 
further supporting the hypothesis that the latter protein 
could have been degenerated as the result of functional 
counter-selection by an active Mod in 5. enterica. 
Despite recent advancements in the structure-function re- 
lationship of Mrr proteins (43), we were unable to 
pinpoint the exact reason for the inactivity of Mrr^^^. 
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Functional incompatibility between Mod™"^ and Mrr'^^"'^^ 

When looking at the phylogenetic distribution in E. coli 
(Figure 6), on the other hand, it seems that homologues of 
mrr ^ and mocf'^'^ are frequently lost or acquired 
throughout evolution. Moreover, in agreement with their 
proposed incompatibility, there is a noticeable anti- 
correlation between both loci, and some E. coli strains 
(e.g. EDIA), in fact, lack an mrr allele while harboring a 
mod allele instead. Importantly, when the mod^ 
gene was cloned and co-expressed it with m;v'^°'*^^ in 
MG1655 mrr, we observed clear evidence for genotoxicity 
(Figure 7). This indicates that, despite bearing only 34% 
amino acid sequence identity with Mod^^^, Mod is 
also capable of rendering the chromosome susceptible for 
cleavage by Mrr'^°'"^ 

Furthermore, similar to LT2, EDIA also became sick 
upon acquisition of pAASlO (carrying the mrr^*^^^^^ 
locus), and readily yielded suppressor mutants (data not 
shown). Interestingly, while half (i.e. five out of 10) of 
these mutants were affected in the mrr^*^'^^^ locus on 
pAASlO, the other half failed to yield an amplicon of 
the mocf'^^^ gene by PCR. In this context, it is note- 
worthy that the mod^^''^ locus is located near a 
prophage integrase gene, perhaps indicating that it forms 
part of a mobile genetic element that can easily be excised 
from the chromosome. 





1600 n 




1 Ann 
14UU - 




1200 - 


'S 


1000 - 








800 - 




600 - 




400 - 




200 - 




0 - 



pACYC184 + 
pFPV25-P5^ + 

pFPV25-Pj;^-modEDiA 



+ 
+ 



-I- 
+ 



Figure 7. Expression of the Mod™'^ MTase triggers Mrr"""*" 
activity in E. coli MG1655. Activation of the SOS response in late 
exponential phase cultures of MG1655 mrr XsfiAwlacZ equipped with 
the indicated combination of plasmids and grown in the presence of 
0.02% arabinose. p-galactosidase activity in the MG1655 XsfiAwlacZ 
background serves as a measure of SOS induction, and is expressed 
as Miller Units. 



Evolutionary antagonism between close homologues of 

Spurred by the clear asymmetric phylogenetic distribution 
of mrr'^'^"'^^ and mod^^^ homologs within E. coli and 
S. enterica (Figure 6), we subsequently examined their 
horizontal transfer index throughout 165 sequenced pro- 
karyotic genomes according to the method of Nakamura 
et al. (42). Although this index only permits detection of 
recent gene transfers from species with distinct gene 
sequence compositions, we nevertheless found that 11% 
(2/18) of Mrr'^°"^^^ homologues and 28% (9/32) Mod^^^ 
homologues were most likely the subject of such recent 
horizontal transfer. Accordingly, we inferred that close 
homologues of mrr^'^'^^^^ and mocF^^ could more gener- 
ally act as mobile elements, and subsequently wondered 
whether the mutual exclusiveiiess of both genes would 
extend beyond strains of E. coli and 5'. enterica. To 
approach this question, we searched 272 completely 
sequenced bacterial and archaeal genomes, and hsted 
those in which Mrr'^°'*^^ and/or Mod^^^ homologues 
were present (Figure 8). Using this fist, the Pearson r cor- 
relation coefficient was calculated for 192 of these 
genomes (with one samphng per genome) and shown to 
be —0.52, which indicated a strong and statistically signifi- 
cant (_P< 0.0001) anti-correlation between both loci. 
As such, the presence of an Mrr homolog tends to be 
commonly paired with the absence of a Mod homolog 
in the same genome, and vice versa. 



DISCUSSION 

During our attempts to rescue the high-pressure induced 
SOS response in S. typhimurium LT2 by introducing 



]yfj.j.MGi655 ^ MG1655, wc obscrvcd that the 
presence of this Type IV REase conferred constitutive 
(i.e. without the need for activation by high pressure) 
genotoxicity in LT2, while a similar toxicity was totally 
absent in MG1655. Since our knowledge about Mrr func- 
tionahty is still limited and fragmentary, we decided to 
dissect the molecular basis of this apparent discrepancy 
between its impact on MG1655 and LT2. As such, we 
were able to show that the observed genotoxicity was 
elicited by the activity of a Type III MTase present in LT2 
(i.e. Mod of the StyLTI system), which most likely makes 
the chromosome susceptible to degradation by Mrr^°'*^^. 
This feature, however, was not limited to Mod^^^ as another 
Type III MTase encoded by E. coli EDIA (i.e. Mod™''^) 
could similarly eHcit Mrr'^*^'*^^ activity. Furthermore, 
because of this functional antagonism, loci encoding 

]y[j.j.MG1655 jyj^ jLT2/EDlA ^^^^ ^y^^^^ ^.^^^jj^ 

counter-select each other. In the same context, we 
revealed a strong and more general anti-correlation 
between close homologues of mrr^'^'^^^^ and mocF^^ in 
the genome database, suggesting that their incompatibility 
might well-extend beyond the genus of Escherichia and 
Salmonella. As such, this study is the first to provide 
evidence for Type III MTases as natural triggers of Mrr 
activity. 



The observed interference between Mrr 



and 



Mod^^^^^^''^ has a number of interesting repercussions. 
As such. Type III MTases only methylate one specific 
strand of its 5-6 bp asymmetric target sequence, indicating 
that Mrr'^"^"'^^^ is, in fact, able to restrict hemi-methylated 
DNA. Moreover, during the writing of this article, in vitro 
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Figure 8. Occurrence and anti-correlation of close homologues of Mrr"^*^'^^^ (blue dots) and Mod^^~ (red dots) in the genome database. 



evidence was published that revealed a similar activity for 
an Mrr-like REase (i.e. MspJl) from Mycobacterium sp 
JLS (44). In addition to hemi-methylation, ah currently 
characterized Type III MTases seem to modify adenine 
(45), indicating that their activity remains undetected by 
other characterized Type IV REases (i.e. McrBC and 
McrA) which specifically target modified cytosines 
(46,47). Finally, our results also indicate that inactivation 
or deletion of the res " gene greatly attenuates Mod 
dependent Mrr'^'^'^^^ toxicity in LT2. In turn, this might 
suggest that the expression and/or activity of Mod ^" 
could be affected by Res^^^, which challenges the 
current assumption that Type 111 MTases function inde- 
pendently of their corresponding REases (48), 

Interestingly, their specificity for methylated or 
modified DNA has allowed the solitary Type IV REases 
to fulfill a number of different activities or roles in the cell. 
As such, it was demonstrated that McrA, McrBC and 
GmrSD are able to restrict invading phages that chemical- 
ly disguise their chromosome with increasingly complex 
base modifications (8,49), Indeed, while McrA and 
McrBC are able to restrict DNA containing hydroxy- 
niethyscytosine (HMC), GmrSD has developed specificity 
for glucosylated HMC, and this differentiation is thought 
to reflect the evolutionary arms race between phage modi- 
fication systems and the host REases, Most recently, 
however, Fukuda et al. (50) demonstrated that the 



deUberate introduction of Type 11 MTases sharing the 
R™C sequence specificity of McrBC could trigger cell 
death in E. coli by eliciting McrBC-dependent chromo- 
some degradation, thereby counter-selecting the establish- 
ment or maintenance of the MTase, As a result, these 
authors postulated a novel function for Type IV REases 
in warding off epigenetic regulation by incoming MTases 
(50,51), The functional and evolutionary antagonism 
between Type IV Mrr and Type 111 Mod systems observed 
in this report seems to corroborate this hypothesis. 
Moreover, aside from being subject to horizontal gene 
transfer, some Type 111 MTases have actually been identi- 
fied as dedicated epigenetic regulators, and have recently 
been shown to direct the expression of a number of genes 
(including virulence genes) in species of Neisseria and 
Haemophilus (48,52,53), 

Interestingly, and fuelling their importance as epigenetic 
regulators, the coding sequence of many Type 111 MTases 
harbors homopolymeric tracts or tandem repeats, 
allowing frequent frame-shifts to occur and making their 
expression prone to phase variation at the translational 
level (48,52), In turn, this results in a phase variable 
ON/OFF switching of Mod-dependent regulons (i,e. 
phasevarions), which is believed to be a strategy that in- 
creases the adaptive potential of populations facing 
fluctuating environmental conditions. Although these phe- 
nomena have not been thoroughly investigated in 
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S. enterica or E. coli, it should be noted that the corres- 
ponding genes contain a number of 6-7 bp homopoly- 
meric tracts. Moreover, we have shown such tracts to 
support phasevariation in the mod^^^ gene upon function- 
al counter-selection of its gene product by Mrr'^°"'^^. 
Furthermore, it is clear from the genome database that a 
number of S. enterica strains naturally contain out-of- 
frame mod alleles. The foil pathogen 5. Gallinarum str. 
287/91, for example, harbors an inactive mod allele 
because of a +1 frame shift in a poly-G tract. Whether 
Mod also contributes to epigenetic regulation in strains of 
E. coli and 5". enterica still remains to be established. 

Aside from the control of phasevarions, however, 
frequent phase variation of Mod activity would have a 
number of interesting consequences with respect to its an- 
tagonism with Mrr. As such, lateral acquisition of Mrr 
could readily select for a population phenotypically re- 
maining in the OFF state by counteracting those cells in 
which the mod allele switches ON. On the other hand, 
transient OFF states might actually have evolved in 
order to allow safe horizontal passage or dissemination 
of mod alleles through strains exhibiting Mrr activity, 
thus mitigating restriction or ceU death. Finally, one 
could even imagine that switching Mod activity in the 
presence of a cognate and functional Mrr would tempor- 
arily allow a small subpopulation to experience bursts of 
DNA damage and increased mutation rates, and to 
behave as transient mutators that can generate genotypic 
variability. The transient nature of increased mutation 
rates is essential, as constitutive mutators eventually 
succumb to the accumulation of deleterious random mu- 
tations (54). Obviously, the vahdity and relevance of the 
possible phenomena stated above require further 
investigation. 

In summary, we present data supporting the functional 
interference and corresponding genomic incompatibihty 
between Type IV Mrr and Type III Mod activity within 
E. coli and S. enterica, and extrapolate that these obser- 
vations might extend to other species. Furthermore, this 
apparent antagonism fits within the recently emerging 
view of Type III MTases as phase variable epigenetic regu- 
lators and the hypothesis that Type IV REases might have 
evolved to counteract the establishment of specific genome 
methylation systems. 
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